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Abstract 

We demonstrate that a highly excited quantum electromagnetic mode strongly 
interacting with a single qubit exhibits several distinct resonances in addition to 
the Bloch-Siegert resonance condition that arises in the interaction with classical 
radiation. The resonance phenomena are associated with non-classical effects: 
Collapse of Rabi oscillations, splitting of the photon wave packet, and field-qubit 
entanglement. The analysis is based on a semiclassical phase-space flow of 2-by-2 
matrices and becomes exact when the photon number tends to infinity. The flow 
equations are solved perturbatively and numerically, showing that the resonant 
detuning of the field and qubit frequencies lie on two branches in the leading 
order, that further split into distinct curves in higher orders. 

1 Introduction 

Physical systems where two-level systems dipole-interact with the quantum electromag- 
netic field are ubiquitous. Once the strong-interaction condition is attained, where the 
qubit-field interaction is appreciably stronger than dissipative environment coupling, 
these systems facilitate the coherent manipulation of qubits. A pioneering setup is 
cavity quantum electrodynamics (QED) [TJ, where atoms interact with electromagnetic 
modes of a high-Q cavity. An analogous process has been realized in superconducting 
circuits, where qubits based on Josephson junctions interact with the electromagnetic 
field in wave guides (2J [3] . 
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The atom-field interaction in cavity QED is often very weak, so that efficient in- 
teraction can only be achieved by making it resonant, in the sense that the frequency 
of the electromagnetic mode is nearly equal to the energy spacing of two atomic levels 
(divided by Planck's constant). Under these conditions the amplitude of atomic levels 
beside the two resonant ones is negligible, and the atom evolves as a qubit; the same 
conditions also imply that the effect of the nonresonant terms in the field-atom inter- 
action Hamiltonian is very weak, motivating the 'rotating wave approximation' (RWA), 
wherein these nonresonant terms are dropped from the Hamiltonian [3]. 

Nevertheless as the interaction becomes stronger, the effects of the non-resonant 
terms becomes noticeable and leads to qualitative changes in the dynamics. These effects 
were first studied by Bloch and Siegert in the context of nuclear magnetic resonance 
where the spin interacts with classical electromagnetic radiation [5]. They showed that 
the resonance curve of the Rabi oscillations of the spin polarization as a function of the 
radiation frequency shifts so that its center is different from the qubit natural frequency; 
the qubit-field resonant detuning is proportional to the square of the interaction strength. 

By the correspondence principle, when the quantum electromagnetic mode is pre- 
pared in a highly excited coherent state the system dynamics should be close to that of a 
qubit in a classical field. This is indeed the case for short times, where the field evolution 
is approximately classical and the qubit Rabi oscillates, but for longer times manifestly 
non-classical behavior takes place, including collapse-revival of the Rabi oscillations [6] , 
and splitting of the electromagnetic field wave packet to create 'Schrodinger cat' states 
|8]. Nevertheless these effects can be understood and described with high accuracy 
using semi classical theories based on flow in generalized phase spaces [91 [TO] . 

Most of the studies of highly excited field cavity QED were carried out in the RWA on 
resonance. When the qubit-field interaction is not very small, the RWA is no longer valid 
and the resonance shifts, and, as shown below, splits into several distinct resonances. 
There is a large literature on the physics of the full Rabi model, see for example [TTI 
[T2l[l3l[lJ[ig[li2[I^ fewer studies have been devoted so far to 

the interaction with highly excited initial wave packets [231 12H [25] . These works used 
semiclassical phase space flow to analyze the collapse and revival of Rabi oscillations, 
wave packet splitting and entanglement, but the semiclassical methods available at the 
time were not sufficient to develop a systematic approximation. 

Here we present a general semiclassical framework for the dynamics of a canonical 
degree of freedom coupled to a two-level system based on the Wigner-Weyl phase space 
representation with matrix- valued symbols [2S] ■ Our theory is based on the propagation 
of matrix-valued phase-space symbols of smooth observables in the Heisenberg picture 
[HI [26], [27] [28] [29]. The symbols are expanded in the semiclassical small parameter e, 
inversely proportional to the mean photon number, and a hierarchy of partial differential 
equations is obtained from the Heisenberg equations of motion. We solve the equations 
to subleading order in e, where the phase-space derivatives are first-order and the equa- 
tions therefore describe a flow in phase space. The leading-order equations describe the 
classical dynamics of the electromagnetic field mode and the Rabi oscillations of the 
qubit in a classical electromagnetic field, while the subleading flow equations describe 
the splitting of the wave packet. 
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The theory developed here is different from that of the well-studied semiclassical 
mechanics of many-component wave function [26, 28, 29J in that the leading term in the 
Hamiltonian is scalar, that is, it acts trivially in the internal degree of freedom subspace. 
The gaps between the bands are then small and the invariant-space decomposition is 
not useful. In a series of papers Bolte and Keppeler developed a semiclassical method 
relying on a flow where the qubit degree of freedom is represented on the sphere [3T1 
|32| [33] . A related method relying on a semiclassical flow in an extended phase space 
has been applied to cavity QED by [MJ [35] . Here we use matrix flow in the standard 
Euclidean phase space to calculate the subleading order in the field dynamics which 
enables the analysis of the splitting of wave packets and generation of Schrodinger cat 
states. The qubit-field back action is studied in the extended phase-space framework in 
a forthcoming work [36J. 

The method presented here is applicable to any initial state where the field is highly 
excited, not necessarily phase-space localized, including thermal states at temperature 
which are high compared with the mode frequency. The phenomena of interest in this 
work, namely the collapse of Rabi oscillations, wave packet splitting, and qubit-field 
entanglement, are captured by the 0(e) equations, leading order qubit and subleading 
field. These effects and the related resonance phenomena take place at times of 0(e^ 1 ^ 2 ) 1 
while the approximation scheme remains consistent for longer times <C s^ 1 . For times 
of 0(e: _1 ) however, the subleading terms become comparable with the leading terms, 
signaling the breakdown of the approximation. Indeed this is the time scale of revival 
of Rabi oscillation that are difficult to capture in semiclassical approximations based on 
phase-space flows [30]. According to [36] the approximation based on the flow equations 
is valid at least for times <C £ -1 ^ 2 , while the results of this paper suggest that this bound 
can be extended to <C e^ 1 , but further investigation of this issue is beyond the present 
scope. 

The semiclassical framework is presented in section |2j and is applied in section [3] 
to the case where the interaction is weak and resonant. The main calculational task 
here is the straightforward problem of the precession of the qubit polarization in a given 
periodic external field. This information is promoted using the semiclassical theory to a 
full description of the collapse of Rabi oscillations, the splitting of the initial field wave 
packet and the generation of qubit-field entanglement. There are natural resonances 
associated with each of these effects. The resonances are all degenerate in the RWA, 
and become distinct once nonresonant terms are taken into account. Our main result is 
the calculation of the resonant detuning for six independent quantum resonance criteria, 
using both weak interaction perturbation theory and numerical solution of the phase- 
space dynamical equations. The resonant detunings are shown in figure [T] as a function 
of interaction strength /i, where the RWA results of zero detuning are recovered when 
fi = 0. The perturbative results are summarized in table [Tj 
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2 Field-qubit interactions 



2.1 Semiclassical Rabi dynamics 

In the dipole approximation, the Hamiltonian of an electromagnetic field mode with an 
annihilation operator a coupled to a qubit with a lowering operator S takes the form 

H = hw\{ota + aa f ) + hu\(S ] S - SS ] ) + %(a + <J){S + 5 f ) (1) 

where a and S satisfy the usual relations [a, a"!"] = [S, S^] + = 1, [a,S] = S 2 = 0, 
and a;, fa/, and % are the mode frequency, qubit energy splitting, and dipole coupling, 
respectively. The Hermitian components of the polarization vector operator S are defined 
by Si = S + , S2 = i(S — S^), and S1S2 = iS%. We do not neglect the nonresonant 
terms aS and a^S^ in the qub it-field interaction. 

We assume that the initial state is a product wave packet |a) <g> |p), where \a) is a 
phase-space localized wave packet, such as a coherent state, with (a| a \a) = a, \a\ ^> 1, 
and |p) is an arbitrary pure qubit state, labeled by the direction of the polarization 
vector p, so that S • p |p) = |p). It is convenient to define a small parameter e such that 
a = (/ \fe with fixed ( and let a — Z/ y/e so that [Z, Z^] = e and 

H = r -f\(Z^Z + ZZ X ) + hv\{&S - S ] S) + h\{Z + Z^){S + 5 f ) (2) 

Since the qubit-field interactions are typically weak, we let A = ^ be finite as e — > 0. 

We can therefore write H = pft+H-S, where ft = \hu(Z^Z+ZZ^), Hi = h\(Z+Z^), 
H 2 = 0, and H 3 = fa/, so that the Heisenberg equations are 

ihd t Z=l[Z,S)] + [Z,H]-S (3) 
hd t S = 2H A S (4) 

Our goal is solve the equations of motion order by order in e. For this purpose we use 
the matrix- valued Wigner-Weyl representation [27], where an operator B is represented 
by a a two-by-two matrix phase-space symbol 6(C) j where ( is a complex coordinate in 
a phase-space with Poisson bracket {(, (*} = —i, such that Z and S are represented by 

z(0 = C , s(C) = CT (5) 

respectively, where a is the vector of Pauli matrices. The product -B1-B2 of two operators 
is then represented by the Moyal product b\j^b 2 = b\b 2 + f {&i, ^2} + 0(e 2 ). The Weyl 
representation of the density matrix of the initial (pure) state is |(1 + <t -p)W a ((), where 
W a (C) is the Wigner function of the field state, so that quantum expectation values are 
given by 

(B) = /^W'(O t r(l±p 6 (C)) (6) 

We assume that the squeezing of \a) remains bounded as e — > 0, so that W is localized 
near ( on the scale y/e. 
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The phase space representation of equations (|3j-|4j) yields equations for the time 
dependent symbols z t and using the Moyal product rule 

hd t z = {z,t)} + e{z,h}-s (7) 
hd t s = 2h A s (8) 



keeping only the leading term in equation (|8j). Assuming small e asymptotics 

z t = zf ) + ez? + ... , s t = sin... (9) 

we obtain 

Hdt4 0) = -idrtf 0) (10) 

ftftsi 0) = 2h<°> A sj 0) (11) 

hdtz? = -id ( *t) {1) - id c *b^ ■ s (0) (12) 

where the asymptotic expansions of the phase-space functions on the right-hand sides 
are defined by 

/ri 0) + ,4 1) + ■■■) = / (0) (C) + ^ (1) (0 + --- 

for / = h, c^*f), and d^*h. Here f) and h are the symbols of Sj and H, respectively. 

The asymptotic approximation ^ is not uniform in time. Nevertheless, the solu- 



tion of equations (10-12) derived below is consistent with the approximation for all 
times <C e' 1 , and captures all the known phenomena of wave packet cavity QED in 
this time interval. We therefore proceed under the conjecture that the validity of the 
approximation can be extended to all orders for t <C e^ 1 . 

Note that unlike the semiclassical equations of motion generated by many-component 



Hamiltonians whose principal symbol is non-scalar J2§J [37] , equation ( 10 -12 ) are not flow 
equations of a classical Hamiltonian. In a companion work [36J we derive Hamiltonian 
equations of motion in an extended phase space that approximate the Heisenberg picture 
dynamics for times of o(e~ l l 2 ). 



2.2 Rabi oscillations and field-qubit entanglement 

The results of the previous subsection mean that CQED expectation values for highly 
excited initial states may be calculated using 
Equation (10) is trivially solved by z t {Q = 
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by solving equations ( |10] - [l2| ) 
so that fj(°) 



-iuit 



hu\(\ 2 and 



(hf'.hf.hf; 



(hfi cos(wt — (/?), 0, \hv) 



(13) 



where \x = A|£| and <p = arg£. Substituting in (11) gives a linear homogeneous system 
for s with periodic coefficients. It is identical with the well-studied equations for the 
Rabi oscillations of a two-level system in a classical radiation field [5j EH]. Here it is 
convenient to apply Floquet theory [39] that implies that the system has three quasi- 
periodic fundamental solutions 



r*(t) 



k = -1, 0, 1 



r fc (t + 



r fc (t) 



(14) 
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The Floquet frequencies are denned modulo u; for concreteness, we choose the branch 
|^fe| < f- Because the RHS of (11) is a cross product, the Rabi frequency Q = Qi = 
is real, Q = 0, and r_ fc (t) = r k (t)* (complex conjugate). 
The solution of (11) can be expressed in the standard way as a linear combination 
of the fundamental solutions (14): Using subscripts a and b for the polarization index, 
the fundamental solutions form a 3-by-3 matrix r a k{t)', denoting the components of the 
inverse of r(t = 0) by r^ 1 we have 

S t ,a = ^W^m, 1 ^ = -b{t>b (15) 
kb b 

defining the orthogonal matrix 0, which describes the Rabi oscillations of the polar- 
ization vector of a qubit in a classical radiation field. The Rabi oscillations are quasi- 
periodic with base frequencies Q and u; f2-periodic oscillations are obtained only in the 
RWA. 

Using ^ we now obtain 

^ = E / ||^(CK nfe W*KtP6 (is) 

kb J 

Recall that r(t) and Q depend parametrically on £ through \i and ip, and this determines 
the ( dependence of Sj. The symbol St(C) therefore varies slowly on the domain of size 
\fe where W a is appreciable, and we may write 

n(c) ~ n(|c|) + (ICI - ICIMICI) = n + (Kl - ICW' , r ak (t; c) ~ r ak (t, () ee r ak (17) 

so that 

(s a > t = £ r afc (t)r fcb V^ I ^W a (C)c^HOT* (18) 
kb J 71 

or, defining the normalized real- variable function w(\(^\ — |£|) = J ^^r^W a (C), 

(sa) t = £ TakWr^w*** J dxw(x)e l ^ tx (19) 



kb 



The integral on the RHS of (19) is the Fourier transform of w. The function w(x) is 
localized near zero on the scale \/e, so its Fourier transform decays on scale t ~ e~ x I 2 . 



Therefore for short times t <C e l l 2 the subleading Q! k is negligible and 

(s) t = 0(t)p (20) 

That is, the polarization of the wave packet Rabi oscillates as a whole with a frequency 
determined by its center. 

For long times t ^> e^ 1 ^ 2 on the other hand, J dxw(x)e iaktx <C 1 so that 

(s) t = Q(*)p, Qab(t) = r a0 (t)r ob 1 (21) 
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and the Rabi frequency f2 disappears from the polarization dynamics. This is the gen- 
eralization of the collapse of Rabi oscillations [6]; unlike the full collapse in the RWA, 
in general the polarization continues to oscillate periodically with frequency u in the 
long-time regime. 

The wavepacket polarization state is an indicator of the entanglement between the 
qubit and the field. Since the field-qubit state is pure, the purity tr Pq Ubit of the field- 
reduced density matrix p qu bit = trgcid p is an entanglement monotone [40J, with tr Pq ubit = 
| when the field and qubit are maximally entangled, and tr Pq Ubit = 1 when the state is 
separable. Since p qu bit = |(1 + ( s ) • CT ) 

tr /4 bit = |(1 + (s) 2 ) (22) 

so that (s) 2 is an entanglement monotone as well. In the short time regime t <C e^ 1 ^ 2 
the wave packet is almost fully polarized since is orthogonal, and the state remains 
almost separable. In contrast, Q is not orthogonal and therefore the field and qubit are 
in general entangled and the degree of entanglement oscillates with frequency u. The 
details of the entanglement dynamics are studied below. 



2.3 Wave packet splitting 



The analysis presented so far was based on the leading order equations of motion (10-11) 
that hold for classical electromagnetic radiation, and the wave packet effects were derived 
by combining them with phase-space integration. We turn next to the subleading field 
flow equation (12). Since d$*t) = (, and d^*hi = A, ^*h 2 ,3 = 0, we have 

dtzjU = -icozi 1] - 2iAs! (23) 



so that has non-trivial qubit action. 



We are interested in the long time t ^> e behavior of z^ , so we use the late-time 
solution (21) for s and Fourier expand r 

oo 

r«fc(*) = Yl h.nae-*^ , t^ = t-~ (24) 



to obtain 

-t 



4 1} = -2iAe- fcrt / dt's^t'y^' 
Jo 



, f _— 1 „ik£lt—i(n—l)u)tm s 



(fe,n)^(0,l) 



All the terms in parenthesis in this expression except the first are bounded. For late 
times, therefore, the first term is dominant and we can write 

Zt = e~ iujt (C + vtn ■ a) (26) 
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up to higher order terms, with 



v = -2i£Ae^f 0il ||r- 1 || , Urol 2 = ^(r^ 1 ) 2 , and n b = r 0fJ 7 1| r 1 1| (27) 

b 

As an operator in the qubit space z t has eigenvalues e~ tuJt ((±vt) with corresponding 
eigenprojections |(l±n-<x). It follows that the expectation value of a classically smooth 
field observable B with symbol b(() is 

(B) = j ^(Ke-^(C + ^))^^ + Ke-^(C-^))^^)^(C) (28) 

which means that the initial wave packet splits into two fragments with statistical 
weights j(l ± p • n) that drift with speeds ±vt in a frame rotating with frequency 

CO. 

The splitting is a necessary condition for appreciable field-qubit entanglement; here 
splitting and entanglement arise on the same time scale t c ~ e^ 1 ^ 2 . The solution (26) is 



valid for time t 3> t c and then vt is much larger than the field uncertainty in the initial 
wave packet and the field state consists of two well-separated sub wave packets. For even 



later times t ~ e^ 1 the subleading term in (26) becomes comparable with the leading 
term signaling a breakdown of the approximation at this time scale. Our analysis is 
therefore limited to the range t ^ s^ 1 . 

2.4 Calculation of the polarization dynamics 

The polarization dynamics Floquet problem is equivalent to an infinite matrix eigen- 
problem for the Floquet function Fourier coefficients [39] 

^ ^ Lna,mb^k,mb ^^k^k,na (29) 
mb 

where the nonzero components of the matrix L are L natna = iun, L n 2, n i = — L n i in2 = u, 
and L„ 3i ( n+1)2 = £( n +i)3,2 = -L n2 ,( n +i)3 = -£(n+i)2,3 = /•*• We present results based on 
perturbative and numerical solutions of the eigenproblem. 

When fj, = the matrix L is block-diagonal. For small /i the Floquet functions and 
Rabi frequency can therefore be calculated by perturbation theory. As our main goal 
is to study resonance phenomena we assume that the qubit-field detuning 5 = v — uj is 
small, and then the spectrum has gaps of order 5. We therefore assume that 5 = O(n) 
and shift a block-diagonal term proportional to 5 to the perturbation matrix, so that 
the unperturbed spectrum consists of all integer multiples of iu, each 3-fold degenerate. 
We calculated the Floquet functions and Rabi frequency to third order in \i using the 
degenerate perturbation theory formulation of [H] . The leading order 0(/x) perturbation 
that includes the interaction between the n = and the n = ±1 blocks is equivalent 
to the RWA. Higher orders go beyond the RWA approximation, and are found to be 
numerically satisfactory up to at least \i = |. 



In addition, we solved equation (29) numerically by truncating the blocks of L with 
\n\ > n max . Very precise solutions were obtained for fi < | with n max = 8 or smaller, 
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Figure 1: The resonant detuning 5 res as a function of the interaction strength /i, both 
expressed in units of the field mode frequency u. The classical Bloch-Siegert detuning 
<5bs is shown in blue, and the quantum detunings 5 tc , 5f c , 5 TC , 5 en , 5 VS , and 5 WS , are 
shown in green, blue, violet, red, green, and violet, respectively. <5f c , <5 VS , and 5 WS , are 
not graphically distinguishable from 8bs, S tc , and 6 rc (respectively), but they are all 
numerically distinct. The graphs of the leading term in the expansion of the resonant 
detunings <5 res = ±\fi 2 are shown in black. The various resonant detunings are defined 
in section [3j 

where it is easy to diagonalize L numerically. Higher values of fi can also be studied in 
this method, but since resonance phenomena are less important when \i is large we did 
not pursue this direction. 

3 Resonance phenomena 

In many experiments the field-qubit interaction energy is much smaller than both the 
qubit energy splitting and the electromagnetic mode energy quantum. In these cases 
the only way to achieve appreciable energy exchange between the field and the qubit is 
to make the interaction resonant, so that the detuning 5 = v — u is comparable with or 
smaller than the interaction frequency scale /i. 

When h/oj — > 0, the rotating wave approximation becomes exact, and the resonance 
condition 5 = implies that the visibility of the Rabi oscillations is maximal, and the 
Rabi frequency Q is minimal for a fixed fi. When fi/u is small but finite the resonance 
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Resonant detuning 


5bs 


5tc 


fife 


<5rc 


<5en 


5vs 


^ws 


Shift coefficient c 


i 

4 


i 

4 


1 
4 


1 
4 


1 

4 


1 
4 


1 

4 



Table 1: Coefficient of the leading order detuning <5 res = c/i 2 for small /x, of the various 
wave packet-qubit resonances, c = ±| in all cases. The resonances are defined in section 

El 



shifts to nonzero detuning 

6 I^ = C (V) 2 + ... (30) 

Bloch and Siegert [5] calculated perturbatively the shift of the minimum of Q as a 
function of £ for a given ^ and showed that c = — \ and f2 res = — Yg(fj) 2 + •••). 
The Bloch-Siegert problem has been revisited several times and high-order terms have 
been calculated [121 S3] . It is straightforward to calculate the Bloch-Siegert shift from a 



numerical solution of the Floquet equation (29); the perturbative and numerical values 



of the Bloch-Siegert detuning are shown and compared in figure [T] and those of Q res are 
shown in figure |2j 

The Rabi frequency is defined by the qubit response to external fields, and the 
Bloch-Siegert resonance therefore characterizes the classical limit e — > 0. We show that 
when e is nonzero, several additional resonances arise from the processes tied to the 
quantum properties of the electromagnetic field. In the rotating wave approximation all 
the resonances occur for zero detuning. In the remainder of this section we study these 
resonances and calculate the resonant detuning of each type. Interestingly, we find that 
while the resonances are all distinct, the leading order detuning of the seven resonances 
(including the Bloch-Siegert resonance) fall under the two cases c = ±|. The values of 
the resonant detunings are summarized in table [T] and figure [T] 

3.1 Collapse of Rabi oscillations 



It follows from equation (19) that the time scale t c of the collapse of Rabi oscillations is 



inversely proportional to d\^\Q = j~d^Q. In the RWA (<9 M fi) = ^— - — is minimal at 
resonance 5 = 0, while for non-zero /x the minimal collapse time, obtained for 5 tc with 



1, is 



tc = t c , RWA {l + ^(-) 2 + ■■■) (31) 



16 \uJ 

The numerically calculated 5 tc and t c are shown in figures [T] and |3j respectively. 

If the initial qubit is a population eigenstate, so that p = e3 is the direction 3 
unit vector, the collapsed population difference in the RWA is (s 3j4 ) = ^ 2+<52 so that 
exact resonance also implies complete collapse of the population. When /i ^ the 
Rabi oscillations fully collapse with (s 3>t ) = identically only for resonant detuning 5 fc 
with c = j. Otherwise the polarization continues to oscillate with frequency u after 
the collapse of Rabi oscillations; still, for detuning 5 rc with c = — \ the polarization 
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0.0 0.1 0.2 0.3 0.4 0.5 

Figure 2: Blue tinted curves: The resonant Rabi frequency f2 res in units of its value 
^rwa = /i in the rotating wave approximation. Red tinted curves: The magnitude of 
the drift speed |i>| of the wave packet fragments in units of Both are shown as a 
function of the normalized interaction strength ^. The perturbative and the numerical 
results are shown in dashed and solid lines, respectively. 
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Figure 3: Blue tinted curves: The resonant collapse time t c in units of the RWA collapse 
time. Red tinted curves: The expectation value of the square of the qubit polarization 
vector s, a decreasing entanglement monotone, averaged over a period 2ir/u at times 
much later than the collapse time t c . Both are shown as a function of the normalized 
interaction strength -. The perturbative and numerical results are shown in dashed and 
solid lines, respectively. 



oscillates with zero mean, (S3 1) 
and S TC are shown in figure [T] 



0. The numerically calculated S{, 



3.2 Field-qubit entanglement 

The resonant Rabi dynamics is characterized in the RWA by maximal entanglement of 
the field and atom, so that (s) = after the collapse of the Rabi oscillations. When 
H is nonzero, (s) does not tend to zero identically for any choice of 5 and continues to 
oscillate with frequency u. We therefore define the entanglement-based resonance where 
the mean square polarization (s t ) 2 = ^ 0+2lT ^ ul dt (s t ) 2 is minimal (see equation (22)). It 
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is obtained for 5 en with c = — \ where 



(s t ) 2 = ^ + --- (32) 

The numerically calculated detuning and mean square polarization are shown in figures 
[T]and[3j respectively. 



3.3 Splitting of a localized wave packet 



The speec 



\v\ of the split wave packet fragments is proportional to ro,i||r || according 

so that \v\ is maximal on resonance. 



-ii 



equation (27). In the RWA ? 01 ||r i — 

The maximum speed is shifted to 5 VS with c = \ for nonzero /x, where 



IC|V 8 + 



(33) 



The numerically calculated S vs and \v\ are shown in figures [I] and |2j respectively. 

For an initial population eigenstate p = e3, it follows from equation (28) that the 
wave packet is split into fragments of equal statistical weight when 113 = 0. In the 
rotating wave approximation 113 = so that the resonance condition implies equal 

weight splitting of the initial wave packet. Equal weight splitting occurs at 5 WS with 
c—\ for nonzero /i. The numerically calculated 5 WS is shown in figure u\ 



4 Conclusions 

The quantum dynamics of the full Rabi dynamics of a qubit coupled to an electro- 
magnetic mode is considerably more complicated than the simplified dynamics in the 
rotating wave approximation, and accordingly less understood. For highly-excited wave 
packets involving a large number of eigenstates, the semiclassical method presented here 
yields accurate quantum expectation values for long times, where the dynamics departs 
radically from that of a qubit in a classical radiation field. It therefore enables us to 
study the collapse of Rabi oscillations, field-qubit entanglement and wave packet split- 
ting effects that characterize the interaction of a qubit and a highly-excited wave packet 
on the basis of phase-space dynamics of two-by-two matrices that operate in the qubit 
Hilbert space. An alternative semiclassical method based on a scalar flow in an extended 
phase space is developed in [3H] . 

The phase-space equations of the full Rabi dynamics are themselves non-trivial. Still, 
they are readily amenable to perturbative analysis near integrable limits and easy to 
solve numerically. The solution suffices for the computation of the quantum observables 
of interest. The focus here has been on the resonant limit, where the interaction between 
the field and qubit is weak and the detuning between the field mode and the qubit 
frequencies is small; the rotating wave approximation becomes exact in the limit that 
both quantities tend to zero together. 
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Bloch and Siegert observed [5] that for a finite-strength interaction with classical 
radiation the center of the resonance is shifted to a finite detuning value. In quantum 
electrodynamics, the picture changes in two ways: first of all, several concepts of reso- 
nance arise, each associated to a certain characteristic observable; secondly, unlike the 
Bloch-Siegert shift, some of the resonances shift to positive detuning. 

We find that the center shift of all resonances take only two possible values in the 
leading order (in qubit-field interaction strength), which are the negatives of each other, 
as summarized in table [T] and shown in figure [TJ This partial lifting of the degeneracy is 
difficult to understand purely on the basis of the phase-space dynamics, but seems more 
natural in view of the approximate transformation [21] of the Rabi Hamiltonian (JT1) to 
RWA form that depends on a single additional parameter. 

Nevertheless, the resonances are all distinct, and split apart once higher order terms 
in the interaction are taken into account, although the splitting of some of the shifts 
is so small that it cannot be observed on the scale of the plots shown here. The com- 
parison between the perturbative and numerical calculations shows that although the 
perturbation theory is ostensibly valid only for weak interactions, the leading order in 
the perturbation series is a numerically excellent approximation for all interaction values 
up to half the field mode frequency, at which point the notion of resonance is no longer 
sharp. 

The main limitation of the present theory is that it is not uniform in time, a property 
common to all dynamical semiclassical approximations. We showed that it is possible to 
extend it consistently for times long enough to cover several non-trivial quantum elec- 
trodynamic effects and analyze the associated resonance phenomena. However, some 
of the wave packet phenomena cannot be captured, notably the revival of Rabi oscilla- 
tions. Such long time effects may be studied using wave packet semiclassical methods 
like those employed in [10J. The advantage of phase-space-flow-based theories like the 
present one is that they apply to any highly-excited initial state and in particular are 
not limited to localized wave packets. 

The experimental observation of the present results is not straightforward, since it 
entails very strong qubit-field interactions as well as long time coherence of excited wave 
packets. On the other hand, the quantum effects arise faster with strong interactions, so 
the necessary coherence time is shorter. A clear indication of a quantum electrodynamic 
resonance would then be a blue shifted resonance curve. Continuous-variable quantum 
information processing in superconducting quantum circuits is a possible application, 
where resonance conditions may arise when seeking resource optimization. 

This work was supported by the German-Israeli Foundation (grant 980-184.14/2007). 
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